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Abstract 

This article highlights interactions of diverse areas: the Heron for- 
mula for the area of a triangle, the Descartes circle equation, and 
right triangles with integer or rational sides. New and old results are 
synthesized. We first exploit elementary observations about circles to 
characterize an arbitrary triangle using three circles. A fourth circle 
brings a certain symmetry - the four radii are exactly the factors of 
the Heron area formula. The four equi-circles also play a role. 

When the triangle is a right triangle, the four circles combine into 
a single tangent cluster. Their centers form a rectangle, the system 
is reflection-congruent to its own dual set, and the dual circles cor- 
respond to the equi-circles (in-circle and ex-circles). The four radii 
satisfy = r\ + r2 + , • r% = r\ ■ r±. Moreover, [—r\,r2, ^3, T3] is 
a quadruple that satisfies the Descartes Circle Equation (DCE): the 
sum of the squares is half the square of the sum. Thus a right triangle 
with legs 2 = ri + r2, 4 = r\ + r% and hypotenuse 2^/h = r2 + r% gives 
this DCE solution: 

[VE-3, VE-l, y/E + 1, VE + 3] 

Pythagorean Triples (PT's) are positive integers that form the legs 
a, b and the hypotenuse c of a right triangle. The basic theory fits 
quite well with the circle method. Convenient radial equation 2r\ 2 = 
(>2 — Ti)(r3 — r\) is equivalent to an old and powerful identity of 
Dickson: 

\{a + b + c) 2 = (c-6)(c-a) 
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Ratios ^r, ^ are the half-angle tangents for the right triangle; the 
reduced numerators and denominators form a Fibonacci-rule sequence 
V. Primitive sequences V and their multiples serve to generate prim- 
itive PT's and their multiples. This includes the ancient standard 
solutions, and new ones. The product of the four numbers of V is the 
area G of the triangle. 

One can replace the in-radius r\ with any ex-radius (7-2, r$, or r^) 
by rearranging and adjusting V . This provides a new and natural con- 
struction of the Barning-Hall ternary tree, which grows all primitive 
PT's from the root triple [3,4,5]. A classical theorem involving the 
nine-point circle provides yet another tree construction. Every primi- 
tive PT when rescaled by factor 1/ G furnishes a Descartes quadruple 
of tangent circles with integral curvatures, and so generates an in- 
tegral Apollonian packing (IAP) containing a rectangle of centers. 
Thus PT's serve to generate an infinite number of in-equivalent inte- 
gral packings. A table of quadruples ends the article. 

1 Going in Circles. 

A plane circle C(A,r) has center A and radius r > 0. In Fig. 1(a) we find 
three circles of arbitrary radii r\,r2,T3 arranged with their centers on a line. 
Each one is (externally) tangent to its neighbors. A final circle is added, 
its center on the same line, tangent as shown (internally), and with radius 
r 4 = r l + r 2 + r 3 . 



A triad of three tangent circles is now constructed, in four different ways. 
First omit one circle, leaving two circles tangent to the third. One of the two 
is then rolled around the third (externally or internally) as in Fig. 1(b) until 
all three are tangent. In Fig. 2 there is one such triad. 




Figure 1: Three circles and their ' sum'. 
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Figure 2: A tangent triad defines a triangle. 



The centers A, B, C are connected to make a triangle, with side a opposite 
vertex A, etc. We easily find (including r 4 for later convenience): 

a = n + r 2 = r 4 - r 3 

b = n + r 3 = r 4 - r 2 (1) 
c = r 2 + r 3 = r 4 - ri 

Can we guarantee that the sum of two sides always is greater than the 
third? That a + b — c for example is positive? Investigating, we are pleased 
to find (2). 

r*i = k(+a + b — c) 
r2 = \{+ a — b + c) 

r 3 = V-a + b + c) (2) 
r 4 = \(+a + b + c) 
= [n + r 2 + r 3 ) 

In other words, the three diameters 2ri, 2r 2 , 2r 3 are precisely those values 
which must be positive in order for a, b, c to form a non-degenerate triangle. 
As a bonus, the value r 4 is seen to be half of the perimeter! 

Before inspecting the remaining triads, we point out in Fig. 2 the in- 
circle of the triangle. It is well known and easy to show that the in-circle is 
orthogonal to all three circles. It is also the circle determined by the contact 
points of the triad. 

Now for the other triads. Oddly enough, joining the centers gives exactly 
the same triangle! Fig. 3 shows that the ex-circles play the same role, mutatis 
mutandis, that the in-circle played in the previous figure. 
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Figure 3: Another way to obtain the same triangle. 



Until now, the fourth circle seemed somewhat unnecessary, but now it 
is seen as essential if the four triads are to form an integrated group. The 
central fact that all four triads give the same triangle is expressed by the 
easy equations (1), (2). 

Let TZ— [ri, r 2 , r 3 , r 4 ] be positive parameters, arbitrary except for r 4 , as 
given. Then TZ makes a rather interesting description of the shape of an 
arbitrary triangle. The triangle area G is neatly expressed as the positive 
square root of the product of all four parameters. 

G 2 = Y[r i = r 1 -r 2 -r 3 -r 4 (3) 

This is a "circular" equivalent to a classical formula, known as Heron's 
formula. In the usual way of expressing this result, we have semi-perimeter 
s and factors s — a, s — b, s — c . However, the "translation" provided by 
(4) is easy to check. 

s = r 4 = |(+a + b + c) 
s-a= r 3 = ^(-a + b + c) , . 

s-b= r 2 = V+a-b + c) ^ ' 

s — c— r 1 = \(+a + b-c) 

Our presentation has the advantage of using "circular reasoning" to give 
each factor on the right-hand-side of (3) a basic geometric meaning. We 
denote by % the tangent triad formed by circles with radii rj,j ^ i. The 
three contact points in the triad determine a "companion" orthogonal circle, 
which is an in-circle 7^ or an ex-circle (7^,7^,7^) for the triangle of centers. 
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Let the circle % have radius denoted by Sj. The in-circle is C(I, S4), and it 
has an in-center I. 

If the center I is connected to triangle corners A, B, C three sub-triangles 
ABI, BCI, CAI result. By adding their areas it is readily found that G = 
r 4 ■ s 4 . Very similar reasoning leads to (5). 

G = n • si = r 2 • s 2 = r 3 • s 3 = r 4 ■ s 4 (5) 

Comparing (5) with (3) we discover that G 2 = f| Sj. But much more is to 
come. 

2 Getting it Right. 

Carbon under great pressure may produce a diamond. In our case, it is 
pressure enough to specialize to right triangles. Let the angle at vertex C be 
a right angle. In other words, the Pythagorean equation (6) holds. 



a 2 + b 2 = c 2 (6) 

Substituting from (1), we find that (6) is transformed to r 2 • r 3 = r\ ■ r 4 , 
hence (7). Also (3) and (5) are simplified to (8) and (9). 

r 4 = 7*1 + r 2 + r 3 , r 2 • r 3 = n ■ r 4 (7) 

G = r 2 • r 3 = n • r 4 (8) 

r a = s 4 , r 2 = s 3 , r 3 = s 2 , r 4 = si (9) 



The identity of the S{ and the r« is striking. It correctly forecasts that the 
case of the right triangle may divulge further special features. The first payoff 
is the coalescing of the four triads into a single diagram. For this we need a 
new section. 

3 Congruent/Dual Tangent Systems. 

First we define the four points 

C = (0, 0), B = (a, 0), A = (0, b), D = (a, b) 
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These form anaxj rectangle, as in Fig.4(a), so that 



a = CB = AD, b = AC = BD, c = AB = CD. 

Here ABC is a right triangle with side a in the X-axis, side b in the Y-axis, 
and right angle C at the origin. 

Triangles ABD, ACD, BCD are congruent copies of the same right trian- 
gle. In Fig. 4(a) below the specific right triangle [a, b, c] = [3, 4, 5] is selected 
for definiteness (and so 1Z = [1, 2, 3, 6]). This diagram can serve for the gen- 
eral case. The four following circles are added (Note: D is not usually on 
circle /C 3 ). 




(a) (b) (c) 



Figure 4: How a rectangle generates four tangent circles. 



/C^CtCU), /C 2 = C(B,r 2 ), /C 3 = C(A,r 3 ), /C 4 = C(D,r 4 ) 

Through (1) we verify easily that circles /Q are pairwise tangent. This 
tangent system will be denoted a. We see that contact points 

l = /L3/L4, 1 2 =/Li/L 3 , 1 3 = /Li/L 2 , ± 4 =/L 2 /L 4 

lie on extended sides AD, AC, CB, BD , and are given by 

Ti = (-r 3 ,6), T 2 = (0,n), T 3 = (n,0), T 4 = (a, -r 2 ) 

In Fig. 4(c) they lie from upper left to lower right on the line A : x+y = r±. 
The remaining two contacts are T 5 = /C1/C4, T 6 = /C 2 /C 3 . The segments 
CT 5 , AT 6 , BT 6 , are radii of lengths r 1 , r 3 , r 2 . They are also hypotenuses 
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of triangles similar to CDB, ABC, and BAC. This makes it easy to calculate 
coordinates. 



—ari — br\ 



ar 3 
c 



br 2 

c 



c c 

We are now in position for a spectacular result. The transformation 
p : (x,y) — > r\ — y, T\ — x is a self-inverse operation that fixes Ti,T 2 , 
T 3 , T 4 and exchanges T 5 with T 6 ! In other words, it is a reflection in the line 

A. Most of the claim is nearly trivial. Showing that p : T 5 < > T 6 amounts 

to verifying the following. 

7*1 b n 
r 2 a + c r 3 



(10) 



Substituting b = r\ + r 3 , a + c = r 4 + r 2 , a = 7*1 + r 2 , b + c = r 4 + r 3 and 
simplifying reduces both to r 2 r 3 = rir 4 , which is just (7). The reflection 
p must therefore take the four-circle system a to another system (3. The 




Figure 5: Reflection Duality of two tangent systems. 

result is pictured in Fig. 5 According to Coxeter [5], the two sets of circles 
are Beecroft duals. It is no great effort to show that this combines the four 
triads into a single system a in a way that exhibits the companion orthogonal 
circles as the dual system (3. 

The four congruent triangles ABC, ABD, ACD, BCD are matched with 
the four /3-circles; the matching circle is an in-circle or an ex-circle by (1). 
Tangency for the in- and ex-circles of the same (scalene) triangle is not pos- 
sible. 
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Theorem 1 The systems a and (3 are Beecroft duals, congruent by reflec- 
tion p : x + y = 7"i, satisfying: 

1. Four of the tangency points lie on the line of reflection. 

2. The other two tangency points are symmetric to this line. 

3. Systems a and /3 share all six tangency points. 
4- Each (3-circle is orthogonal to three a-circles. 

There is something very satisfying about the eight interlocking circles 
illustrated in Fig. 5. The complex of relationships here relies on the very 
elementary equations (1) and (7). The alert reader may be curious at this 
point how those equations are related to the celebrated identity of Descartes 
(sometimes called the Soddy formula), which involves the sizes of four mu- 
tually tangent circles. To this identity we now turn. 

4 The Descartes Circle Theorem. 

We rely heavily on [15]. Stated in essence in a letter (1643) by Descartes, 
this theorem was put more broadly and proved by Steiner (1826) and then 
others. Nobel Physicist F. Soddy made a delightful poem, and obtained a 
generalization to n dimensions (1936). 

First we define the oriented curvatures ki which correspond to our radii 

n. 

K=[k u k 2 , k 3 , h} = {± ± i -i] 

Basically the curvature is the inverse of the oriented radius. The oriented 
radius is either r or — r, where the unsigned radius is r > 0. In our context 
the negative sign is used exclusively for the fourth circle, and that is because 
it encloses the others. 

The radial equations (7) can be rewritten as the following analogous 
equations (11). 

^2^3 + ^1^4 = 0, k 4 + k 1 = k 2 + k 3 (11) 

The Descartes circle equation (DCE) is sometimes attributed to Soddy. It 
states that the sum of the squares of the oriented curvatures is half the square 
of the sum (12). 
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[k\ + kj + k 2 3 + kl\ = l[kt + k 2 + k 3 + h} 2 (12) 

This equation follows fairly easily (details omitted) from (11). The burden of 
the extra variables ki can be dispensed with. Since (12) is homogeneous, we 
can use (8) to replace /C by G7C = [r 4 , r 3 , r 2 , — rj. Equally, the triple/triangle 
[a, b, c] is rescaled to -^[a, b, c] to replace K by G/C. With given triple ±[3, 4, 5] 
we have 

ft =±[1,2,3,6], /C= [6,3,2,-1] 
This works. So does /C = [20, 5, 12, -3] from ±[15, 8, 17]: 

(36 + 9 + 4 + 1) = ±(6 + 3 + 2 - l) 2 

(400 + 25 + 144 + 9) = ±(20 + 5 + 12 + 3) 2 

Exercise: Find the four right triangles that are the sources for each of these 
solutions of the DCE. 

[\/2 + 2, V2,\/2,\/2-2] [V3 + 3,V3 + l,V3-3, V3- 1], 

[V5 + 3,V5 + l,V5-l,V5-3], [15,10,3,-2] 

We will return to the topic of Descartes circles at the end. There is one other 
theme we need to pursue first. 

5 Pythagorean Triangles. 

One might happen to notice that many of our examples are right triangles, all 
of whose sides are whole numbers, like [3,4,5] or [5, 12, 13]. A triple [a, b, c] 
of positive integers satisfying (6) is called a Pythagorean triple. 

Our equations (1), (2) show that when the radii are integers, so are the 
sides of the (arbitrary) triangle. They also imply that when the sides a, b, c 
are integers, and at least one radius is an integer, then the other three radii 
are the sum/difference of two integers. 

It remains to consider when the sides are integers, but the radii are half- 
integers, by (2). For instance, the 1-1-1-equilateral triangle has 

TZ = [± ± ± 21 

/v 1-2' 2' 2' 2-1' 

Doubling the sides and the radii then gives an integer solution in which every 
radius is odd. 
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For right triangles this cannot occur by (7). Odd numbers are of two 
types: either Ak + 1 or 4k — 1. How many rj are of the second type? The first 
equation of (7) has a "balance" that requires an odd count, but the second 
equation requires an even count. Thus we conclude that any Pythagorean 
triple [a, b, c] must correspond to integer radii 7Z and conversely. 

Besides, if the radii are multiplied by a factor k > 1, the same factor is 
applied to the sides, and vice versa. This means that we can stick to primitive 
triples, which have GCD equal to one. Clearly this is true for the triple if it 
is true for 7Z, and conversely. 

However, our desire to employ 7Z as the ideal set of parameters for a 
triangle is severely challenged in the right triangle case, since (7) must 
be satisfied! One strategy that can be used is to examine the restricting 
equations very closely. From (7) we remove the redundant r^ and study 
fi ■ (ri + r2 + r%) = r2 ■ r%. The first step is to multiply out and rearrange the 
terms: r 2 • r^ — r\ ■ r2 — r\ ■ r^ — r\ 2 = 0. This very nearly factors (the last sign 
is wrong), but completing the product anyway gives us a useful equation. 

2n 2 = (r 2 - ri)(r 2 - r{) 

or 

\{a + b - cf = (c - 6)(c - a) (13) 

The second equation is obtained by substitution from (2). If we introduce 
variables e,m,n for the three parentheses of the second equation, then (14) 
results. 

a = m + e, b = n + e, c = m + n + e, e 2 = 2mn (14) 

The essence of (14) was introduced a long time ago, by L. E. Dickson 
[6]. We agree with Gerstein [9] that it has been long overlooked; witness 
its rediscovery by Teigen and Hadwin [20] and others. We deal with this 
equation more fully elsewhere [3], showing that it can function as the source 
(almost effortlessly) of nearly all good and basic facts about Pythagorean 
triples. But perhaps the same is true of the circles approach in this article. 

Klostergaard [13] and McCullough [17] use equivalent variables and equa- 
tions to carry out enumerations. Lonnemo [16] uses it to construct the cel- 
ebrated "family tree" (introduced in 1963 by Barning [2] or in 1970 by Hall 
[11]). Preau [18] uses a similar transformation to rediscover the tree. Alperin 
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[1] also rediscovers it. We shall shortly give entirely new methods employing 
circles! 

First we sketch a general algorithm for selecting the radii for a primitive 
triple. One example will suffice. Let n = 3, say. Write 2ri 2 as a product of 
prime powers and one: 18 = (1)(2)(9). With (1) and (2) in different camps, 
we may assign (9) to either camp; thus 18 = (1) (18), or 18 = (9)(2). The 
original 'three' is then added to both factors, and the sum r 4 is computed. 
Thus r\ = 3 leads to two possibilities: 

TZ = [3, 1 + 3, 18 + 3, 28] -> T = [ 7, 24, 25] 

TZ = [3, 9 + 3, 2 + 3, 20] -> T = [15, 8, 17] 

If we eliminate the restriction to primitive triples, it is easier to start by 
picking any value for 

m = c — b = 7*1 — r 2 . 

Then write 2m = f 2 g, where f 2 is the largest possible square factor (and 
consequently g is square-free). Finally, n = gh 2 , and e = fgh. This is a small 
simplification of McCullough's enumeration technique [17]. Did we mention 
that e = 2r*i is a diameter? (Yes, in the first section.) 

We adopt the convention that when [a, b, c] is a primitive triple, side a 
is odd. Equivalently, the radii alternate between even and odd: we already 
know that at least one radius is even, but the second equation of (7) guaran- 
tees a second, and the first equation tells us two are odd (three even implies 
four even, contradicting primitivity). The second equation (again) makes one 
of 7*1,7*4 even and one odd, likewise for r 2 ,r 3 . We exchange r 2 ,r 3 if needed 
(also a, b) to get alternation between even/odd. 



6 Raising the Standard. 



As further proof of the utility of circles, we now extract the age-old standard 
solutions for Pythagorean triples. The key is a "minor" lemma from Section 
3, namely (10). Here it is again, enlarged a bit, courtesy of the identities, 
(6), (7). 

H r$_ b c~a 

T2 ta a+c b 

(15) 

ri V2_ a c—b 

T3 T4 b+c a 
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We avail ourselves of the opportunity to limit consideration to primitive 
triples. Hence the sequence (rj) has no non-trivial common factor, and alter- 
nates in parity. Consequently, the second group of equal fractions includes 
one with odd numerator and odd denominator. 

We denote the reduced fraction in the first group by q/p, and the reduced 
fraction in the second group by q'/p'. The pair (q,p) and the pair (q',p') are 
both relatively prime, and the last is a pair of odd numbers. 

It then works out that r\ = xq, r 2 = xp, and r 3 = yq, r 4 = yp, where x 
is the GCD of r±, r 2 and y is the GCD of r 3 , r 4 . We note that any positive 
factor k of x, y is a factor of all four radii. Hence the pair (x, y) is relatively 
prime. Finally ^ = ^ = | = ^ establishing x — q', y — p' and hence (16). 

v x = qq', r 2 = pq', r 3 = qp', r 4 = pp' (16) 

The key numbers q, p, q', p' are in fact pairwise relatively prime, since any 
factor k that divides two of them must divide three radii, and by (7) all four. 
There is much more. First of all, a right triangle with legs c + a, b is easily 
shown to have one acute angle equal to half an acute angle of right triangle 
[a, b, c]. Checking details shows that q/p and q'/p' are the half-angle tangents 
for [a, b, c}\ 

If we substitute from (16) into equation t\ + r 2 + r 3 = r 4 the result can 
be "solved" in either of the following two ways. 

q _ p'-g' i_ _ p-q 

p p'+q' ' p' p+q 

Since we know that q', p' are odd, we can put qz = \{p' — q'), pz = \{p' + q')- 
But then z(p — q) = q', z(p + q) = p', implying that z = 1. And so 

p-q = q', p + q=p', - q') = q, \{p' + q') = p. (17) 

We can even conclude that in the pair (q, p) one is odd, and one is even. 

Using (17), it is amusing to present the sequence V = [q',q,p,p'] as a 
Generalized Fibonacci sequence! In other words, p = q' + q and p' = q + p. 
One may pick, in fact, q' to be any odd positive integer, and q any positive 
integer which is relatively prime to q'. For example: 

V = [1, q, 1 + q, 1 + 2q] -> K = [q, 1 + q, q(l + 2q), (1 + q)(l + 2q)\ 

^T= [l + 2g, 2g(l + g), 1 + 2q{\ + q)] 
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or with q not divisible by three, 



K = [3, q, 3 + g, 3 + 2g] = [3q, q{3 + q), 3(3 + 2g), (3 + g)(3 + 2g)] 

-> T = [3(3 + 2g), 2g(3 + g), 9 + 6g + 2g 2 ] 

More generally, we get, combining [a, 6, c] = [7*1 + r 2 , 7*1 + r 3 , r 2 + r 3 ] with 
(16), (17): 

[a,b,c\ = [p 2 -q 2 , 2pq, p 2 + q 2 } (18) 

[6,a,c] = [I(p' 2 + g' 2 ), p'q\ \{p' 2 + q' 2 )] (19) 

[a, 6, c] = [pV , 2pg, pg' + qp] (20) 

[a, 6, c] = [pV, 2pg, pp' - qq'\ (21) 

Here (18), (19) are the two standard solutions, the first in terms of q,p 
alone, the latter (with a,b switched) in terms of q',p' alone. We dare to 
add two mixed forms (20), (21), which seem to us to be the most generally 
useful! (Also, note the area is simply G— qq'pp')\ 

Both standard solutions are found in Sierpinski [19], in Dickson [7], in 
Eckert [8], and doubtless other places. One or the other, usually with an 
added factor of k > to allow for non-primitive solutions, are found in 
sources too numerous to list. The mixed forms are close in form and spirit 
to "Generalized Fibonacci" sequences, for example Koshy [14] or Horodam 
[12], but are bit more symmetric than any that we have seen. The results of 
this section are recapitulated in the following theorem. 

Theorem 2 (Radius Decomposition). Let [a,b,c] be a primitive triple 
with radii TZ = [r 1 , r 2 , r 3 , r 4 . ] Then parameters p, q, p', q' are determined 
such that, (16), (17) are valid, and also (18), (19), (20), (21). Moreover, the 
fractions of (15) are the half-angle tangents of the right triangle. 

7 Tree Rings (More Circles). 

Recall that V = [q r , q, ■ ■ ■] can be completed by "Fibonacci addition", V = 
[q', q, p, p'], and we only require that q' is odd, and q',p are relatively 
prime. The four radii r 1; r 2 , r r 4 are found by multiplying either end value 
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by either middle value- this is exactly the result of expanding the product 
(q' +p'){q +p) Consider the possibility of modifying V by exchanging q with 
p, or q' with p', or both. The result is then one of these. 

Vi = W,Pr--l T>2 = W,q,---], V 3 = [p',p,---}. 

The final two values are replaced by computing with Fibonacci addition, 
and in each valid sequence is obtained. 

We ask, what is the in-circle radius for Vi, V2, V{1 It is found as the 
product of the first two members, namely q'p, p'q, p'p. But these three values 
are just ^3, (the ex-circle radii from the starting case V)\ In effect, we 
produce V\, V2, or V3 from V by "promoting" one ex-circle to the position 
of in-circle! Such promotion always increases the size of the in-circle. Take 
the simplest example V — [1,1,2, 3] which has radii TZ = [1, 2, 3, 6]. Promote 
r 2 to get V\ = [1, 2, • • •], or promote r 3 to get V2 = [3, 1, ■ • ■], or promote 
r 4 to get V3 = [3, 2, • • •]. Continuing with V3, we find V3 = [3, 2, 5, 7], then 
multiply out (3 + 7)(2 + 5) = (6 + 15 + 14 + 35). Finally, add (6 + 15), (6 + 14), 
(15 + 14) to get the triple [21, 20, 29]. 

The cancellation of a "promotion" is a "demotion". We claim that any 
V other than V — [1, 1, 2, 3] has a unique demotion! Demotion decreases the 
in-radius. 

Theorem 3 (Demotion) . Let V = [x,y, •••] be a valid sequence other 
than [1,1,2,3]. Then there are three uniquely determined Fibonacci-rule se- 
quences of form [x, *,y,* ], [*, y,*,x], [*, *,y,x]. Exactly one of these has 
all entries positive. This is the demotion V~ of V. Thus V can be demoted 
in exactly one way. 

The proof of this theorem is so direct, we need only sketch the demotion 
algorithm and leave the rest to the reader. If x = y then V = [x, x, 2x, 3x] 
= [1, 1, 2, 3], so assume x < y or x > y. 

If x < y then [x, *, y, *} = [x, y — x, y, 2y — x] works and the other two 
do not. If x > y take [*,*,y,x] and "back up" to [*,z,y,x\. Note z = y 
is excluded, since x = z + y is odd! Then, if z > y, reorder to [*,y,z,x] 
(sort into ascending order). Finish backing up. The result is then either 
[*, *> 2/> A = i 2 V ~x,x-y,x, y] or [*, y, *, x) = [x - 2y, y,x-y,x]. Note the 
initial element is always odd. 

Unlimited demotion (infinite descent) is impossible, for the in-radius is 
always decreasing. Therefore, any triple can be obtained in one way by a 
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sequence of promotions, starting with T = [3,4,5] , that is, V = [1,1,2,3]. 
The set of primitive triples then has the structure of an infinite ternary tree. 
At the top is the "root" % = [3,4,5]. Each "node" is a primitive triple with 
three immediate successors. These successors are (A) left, (B) middle, or 
(C) right respectively as the promoted circle has size r^,r^r2- For the root, 
this means [15, 8, 17] is left, [21, 20, 29] is middle, and [5, 12, 13] is right. The 
repertoire of simple connection equations in Theorem 2 makes it easy to 
go from one to another. For illustration, we compute P-sequence [3, 2, 5, 7] 
and 7£-sequence [6, 15, 14, 35] for the middle triple T = [21, 20, 29]. Readers 
should produce enough examples to demonstrate this for themselves. 

Corollary 1 (B arning/ Hall Family Tree) . The set of primitive Pythag- 
orean triples can be seen as the nodes of a ternary tree. Demotion moves up 
one level to the "parent" and promotion moves down one level, producing 
three immediate successors "children". Only one triple [3, 4, 5] has no parent. 

A certain amount of tedious computation is all that stands in the way 
of expressing promotion/ demotion as a transformation of triples [a, b, c], or 
a transformation of half-angle tangents (q,p) , (q',p')- See for example [3]. 
Here the following will suffice. 

Corollary 2 (Half-angle tangents) . Let | be either half-angle tangent 
for [a, b, c], but not ~ or |. The fraction is demoted to ±(^^) ± . The first ± 
is chosen to avoid negative quantities; the second ± to avoid fractions greater 
than one. Examples: 

_3 3 3 3 2 5 5 5 4 

10 4'8 2 3 '6 2 4 5 

Corollary 3 (Reversed Barning-Hall Transformation) . If triple [a, b,c] 
is not equal to [3,4,5], it may be demoted to 

\a-2e\, \b-2e\, c - 2e = | - a - 2b + 2c|, \-2a-b + 2c\, -2a - 2b + 3c 

Example : [33, 56, 65] -> [ |33 - 48|, 56 - 48, 65 - 48] = [| - 15|, 8, 17] . 

Comments: See (14) for e = a + b — c. As Lonnemo [16] points out, (14) 
works just as well with e changed to — e. This clue leads to the transformation 
of Corollary 2 except for the absolute value operations. 
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The reader may wish to work out the corresponding promotions implied 
by these corollaries: ratio | promotes to 2p+q' p^W' triple K ^ c 1 

promotes to 

[±a±2b+2c, ±2a±b + 2c, ±2a ± 2b + 3c] (22) 

where in all three terms the pattern is the same one: ( — h +) or (+ + +) or 
(+"+)• 

If the triple is taken as a column vector X, then promotion is just a matrix 
linear transformation: X — > AX or BX or CX. Here B is square with (1, 1, 3) 
on the diagonal and 2 elsewhere, and changing the signs in the first (the 
second) column gives A (gives C). 

We have already cited the independent discovery of the ternary tree, seven 
years apart, by Barning and Hall. Eckert [8] alludes to the ideas of Corollary 
2. See Dickson [7] for the fact that Fermat essentially found transformation 
V — > "P 3 , but (apparently) not either of the other two. Of all the different 
ways to construct the tree, we believe our circle method is the simplest, and 
the most transparent. 

We have yet one more original way to construct the tree via circles! The 
key is a beautiful circle from classic geometry, studied by Feuerbach and 
many others. We think of it in this context as the "family circle." 

In preparation, we recall Fig. 4 and Fig. 5. The four tangent circles in 
cluster p are "attached" as equi-circles to congruent right triangles ABC, 
ABD, CDA, CDB. Focus on ABC. The equi-circles for ABC are C(Ii,Ti) 
where 

Ii = (ri,ri), I 2 = (r 2 ,-r 2 ), I 3 = (-r 3 ,r 3 ), I 4 = (r 4 ,r 4 ). 

The four symmetries of rectangle ADBC are obtained from the following: 

a : x — > x, r : x —> a — x 

V -> b-y, y -> y 

These are reflections, vertical and horizontal. Recall also 

P ■ (x,y) -> (ri - y, ri -x). 

The four a-circles K = C(C, ri), C(B,r 2 ), C(A,r 3 ), C(D,r 4 ) are mapped to 
the above equi-circles by K — > pJC, rpJC, crpK., arpK. 
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The center M = (|, |) of the rectangle is also the midpoint of diagonals 
AB and CD. Circle C(M, |) is seen to be the circum-circle of the rectangle, 
and also of the four right triangles. Let F = y = (|, ~). Circle C(F, |) has 
CM (one radius of the circum-circle) for a diameter. It is also a circum-circle 
for a "shrunken" rectangle A'B'C'D': We note that A'B'D' are the midpoints 
of the sides of ABC. Hence this last circle with center F is the nine-point 
circle M studied by Feuerbach and others. 

In Fig. 6 the equi-circles, the circum-circle, and the nine-point circle are 
shown. Note point L, the foot of the altitude from C. The AB diagonal line, 
and the altitude CL have equations bx + ay = ab, ax = by. Hence the foot 
of the altitude is 

t / ab 2 a?b \ 

L ~ I c 2 ' c 2 I' 

Origin C is the foot of the other two altitudes, and also a point of tangency 




Figure 6: In-circle, ex-circles: all tangent to nine-point circle N '. 

between the circum-circle and its "half" brother M . As diameter CM of M 
is also hypotenuse for right triangle CLM, the point L is on circle M. It is an 
interesting exercise to verify this by finding the side lengths, and using the 
Pythagorean theorem. A helpful datum for this pursuit is the Pythagorean 
right triangle with hypotenuse c 3 and legs equal to the absolute value of the 
real and imaginary parts of the third power (a+bi) 3 . In Fig. 6 circle M seems 
tangent to the four equi-circles - no surprise to geometers! We shall verify 
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this beautiful result of Feuerbach, and in the process, make a wonderful and 
serendipitous discovery! The key is to compute the (vector) difference of the 
centers of two circles, and the sum/difference of the radii. Should the values 
(aside from sign) turn out to be the legs and hypotenuse of a right triangle, 
the circles are tangent. Thus 

I4 — F = (r 4 - f , r 4 - J) = J(a + 2b + 2c, 2a + b + 2c) 

Sum of radii: r 4 + f = ±(2a + 2b + 3c). 

The "fingerprints" of (22) have shown up (middle case) — the segment I4F 
is the hypotenuse of a quarter-sized copy of the middle child on the tree! 
By similar reasoning, the segments I2F, I3F go with the other two children 
(quarter-sized). 

Fig. 7 magnificently illustrates the production of the three "infants". But 
it does nothing to explain why the ex-circle used to construct one of the 
"infant" triangles is exactly four times the size of the in-circle of that triangle, 
a fact we know because of the study of promotion at the beginning of this 
section. 



v P.4,5] 




Figure 7: Children of [3, 4, 5] emerge quarter-sized from tangencies. 

But what about the other two circles tangent to circle A/", namely the 
in-circle with center I and the circum-circle with center M? Segment FM is 
the hypotenuse of a quarter-sized copy of the main triangle, and FI is the 
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hypotenuse of the "parent" of the main triangle (all quarter-sized, except, of 
course, the root triple [3,4,5], which has a degenerate parent [0, 1, 1]. Truly 
the nine-point circle of a Pythagorean triangle is a "family" circle! 

One us (Price) was studying the results of Fig. 7 and found a convenient 
way to produce the immediate family (one predecessor and three successors) 
of a given triple T = [a, b, c]. Compute t = 2 (a + b + c) and then check that 
[t — a, t — b, t + c] is a positive triple with a larger hypotenuse, with the same 
GCD (so both are primitive if either is). Note that t is double the perimeter 
of triangle T . 

Now suppose that exactly one of a, b, c is negative. The three cases may 

be distinguished by the sequence of signs: ( — h +), (H h), (+ H — ). We 

add (+ + +) for the all-positive example of the previous paragraph. For 

( — h +), (H h), t/2 is not the perimeter of course, but the result is still a 

triple with larger hypotenuse. In the last case (+ H — ) we must take absolute 
values to get a positive triple: \t — a\, \t — b\, \t + c\. This time, the hypotenuse 
is smaller. 

Changing two or three signs of a, b, c only duplicates the four cases already 
seen, plus some added sign changes. We can now recognize the left, middle, 
and right child of a triple as the same as the outcome of cases ( — h +), (+ + 

+), (H h), and also recognize the fourth case (++-) as the predecessor! In 

hindsight, the method is close to the generalized Fermat method of Eckert 
[8]. 

Finally we are ready to wrap up the fireworks with a finale, one with lots 
and lots of circles! 

8 A Descartes-Pythagoras Reunion. 

In Section 4 we derived the Descartes Circle Equation for a very special 
case. With [15] as reference, we now claim it for any set of four circles, 
tangent at six distinct points. 

Our immediate goal is to examine certain integer solutions. We already 
have the basic solution K,' = [r^,r 3 ,r 2 , —ri] belonging to rescaled triangle 
-g[a,b, c]. Provided that [a, 6, c] is a primitive Pythagorean triplet, integer 
quadruple K! is also primitve. 

The triad of tangent circles with curvatures r^r%^r2 are tangent to the 
surrounding circle with curvature — r 1; but are also all tangent to a different 
circle placed in the midst of them. Let that smallest circle have radius 
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G/r (or radius r for [a, b, c]). The curvature (positive) turns out to be 
G/r = Gk = 4r 4 - r x . For the triple [3, 4, 5] we find Gk = 4(6) - 1 = 23. 
It is easy to check that [6, 3, 2, 23] satisfies the DCE. 

Let f,g, h be any three curvatures, /, g positive, and y + ^ + ^>0. A 
corresponding triad of three circles tangent at distinct points may be con- 
structed; the third circle encloses the others in case h is negative. We may 
substitute /, g, h into (12) to obtain a quadratic, x 2 — Sx + T = 0, with two 
real solutions xi, X2, such that S = x\ + x% and T = x\ ■ xi- 




Figure 8: A different "in" circle for the triangle |[3, 4, 5]. 

It is both clear and commonly known that two different fourth circles 
exist, tangent to the given three. Evidently X\ y Xi are these two curvatures. 
Exchanging one circle for the other is called Descartes reflection. We also 
find 

Xl + x 2 = S = 2(f + g + h). (23) 

Thus [6, 3, 2, -1], [6, 3, 2, 23] satisfy (23), since 23 + (-1) = 2(6 + 3 + 2). 
We see from (23) that Descartes reflection can be applied to any one circle 
out of a set of four tangent circles. Moreover, if all four circles have integer 
curvature, so does the new circle produced by reflection! 

Descartes reflection is a special case of circular inversion. Consider the 
in-circle (dashed) in Fig. 8. Reflection in this circle must fix the three circles 
to which it is orthogonal (2,3,6), and swap the other two (23, —1) because 
tangencies are preserved. Much more information is available in [15] or [10], 
including simple means of locating the centers of reflected circles. Our inter- 
est in Descartes reflection is more limited in this paper. 
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Suppose we start with Fig. 8 rotated 90° clockwise. Then by line reflec- 
tions, or by D-reflections of individual circles, we get Fig. 9(a), an arrange- 
ment of nine circles. 




Figure 9: Descartes quadruple [6, 3, 2, —1] leads to a packing of circles. 

Fig. 9(b) is an enlarged view of the upper left of 9(a). When circles fill 
the gaps (a, b, c, d); the curvatures 11, 14, 23, 2 are found by using (23). 

3 + x = 2(6 + 2-l), 2 + x = 2(6 + 3-1), 
-l + x = 2(6 + 3 + 2), 6 + x = 2(3 + 2-1). 

Since all four calculations start from [6,3,2,-1], it is possible to package 
them together. Replicate twice the sum [20,20,20,20], and also multiply 
each curvature by three, to get [18, 9, 6, —3]. Subtract as vectors: 

[20 - 18, 20 - 9, 20 - 6, 20 + 3] = [2, 11, 14, 23]. 

Each of these values can replace the original value in the same position. 

[2,3,2,-1], [6,11,2,-1], [6,3, 14,-1], [6,3,2,23] 

Clearly this process can be continued indefinitely, filling the interior of the 
"negative" circle with an infinite "foam" of circles of increasing integer cur- 
vature. This is an integral Apollonian packing (IAP). In both [10] and [15] 
is a beautiful picture which carries our Fig. 9 to many more levels. We will 
shortly present a detailed picture of our own. 

The authors of [10], [15] offer no connection between integer packings and 
integer right triangles. But in Fig. 10(a) see how the centers of the largest 
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Figure 10: Packings from Pythagorean rectangle, or symmetric starting 
point. 

circles form four rectangles. The shaded right triangle, similar to [3,4,5], is 
for us a nucleus for the packing. 

There is little that is privileged about the |[3,4,5] triangle. We can get 
different Apollonian packings from any rescaled primitive triple ^[a, b, c]. It 
will always be an IAP with two rectangles, and bilateral symmetry. The 
children of [3,4,5], rescaled, are 

£[15,8,17], ^[21, 20, 29], £[5,12,13]. 

This gives us three more DCE solutions K,' in which we now sort the 
entries to put it in the form [—a, b, c, d] with 0<a<6<c<<i: 

[-3,5,12,20], [-6,14,15,35], [-2,3,10,15]. 

The problem is that different sequences X7 may belong to the same IAP. 
Given two sequences, each determines a unique IAP, but how can one tell 
if the packings are the same or distinct? A good question, but first we show 
how to generate many more packings! That will bring us near to the end of 
our study. 

As we prepare to "pack it in" we return "full circle" to a figure reminiscent 
of Fig.l, namely Fig. 10(b). Suppose that r, s,t are positive integers, and 
that three tangent circles with oriented curvatures r, s, —t have their centers 
on a line. Since the radii form a sum, the (absolute) curvatures must make 
a harmonic sum: 

1 = 1 + 1 

t r ~ s 

We desire that the Descartes equation applied to KJ = [—t, r, s, x] will have 
two identical integer roots x : x = X\ = x^. In that way we can generate an 
IAP with bilateral symmetry. The harmonic sum can be written in the form 
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(r — t)(s — t) = t 2 , which resembles (14), and has this Diophantine solution 
(easily obtained). 

r = km(m + n), r = kn(m + n), t = kmn. 

Here k is the GCD for r,s,t, assuming that m,n are relatively prime. 
We want primitive solutions, so we put k = 1. To avoid duplication by 
symmetry, we put m < n (let m = n = 1 be designated as isolated). Again, 
serendipity! The condition for identical roots x = u in the Descartes equation 
is automatically fulfilled, so that (23) becomes u = r + s — t, thence (24). 

r = m(m + n), s = n(m + n), t = mn, u = m 2 + mn + n 2 (24) 

Fig. 9(a) clearly shows that the isolated solution /C = [—1,2,2,3] and 
solution KJ = [—1,2,3,6] are in the same packing (6 reflects to 2 and vice 
versa). We therefore consider m < n. Suppose we have arbitrary primitive 
triple [a, b, c] = ri+r 2 , ri+r 3 , r 2 + r 3 . Then generator pairs (q,p) and (q',p') 
can both take turns as m,n in (24), in which case the resulting Descartes 
solutions can be written as (25). 

K> = [f, r 3 , r 4 , £±6], K! = [-a, 2r 2 , 2r 4 , 2c + a] (25) 

Take triple [7,24,25] for instance. From the generators — = |,y we derive 
by (24) or (25) two bilateral solutions [12,21,28,37] and [-7,8,56,57=]. 
Here k == indicates that a curvature D-reflects to itself. These agree with 
(25). We also have the more natural KJ = [—3,4,21,28]. 

There is another type of bilateral symmetry to consider. It splits via 
parity into two subcases, illustrated by the DCE solutions [—4, 8, 9, 9] and 
[—3, 5, 8, 8]. A fairly direct analysis, similar to the above yields the following 
algebraic solution for [—a, b, c, c]. 

a = 2k(t-k), b = 2k(t = k), c = t 2 

In this solution, we suppose that | > 2 is a reduced fraction. If t is even, 
then divide a, b, c by two. 

If the packing is identified by the simplest quadruple it contains (largest 
possible circles), called the root, then quadruple [—3,4,21,28] given above is 
not a root. Producing the root for this case is easy. Just reflect the largest 
number three times. 

[-3,4,21,28] -> [-3,4,16,21] -> [-3,4,13,16] -> [-3,4,12,13] 
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More generally, let /C = [—it, r,s,t\. Let —u + r + s = v. Then if t exceeds 
v, D-reflection replaces t with a smaller integer. Repeat until the the process 
ends, either with [—a, b, c, d == ] in which d reflects to d, or with [—a, b, c, d ><e ] 
where d < e and d reflects to e. 

In [10] it is shown that this reduction always leads to a root quadruple 
K, defining the IAP. Our process finds many - see the table at the end of 
the article. One of the simplest packings (see Fig. 11) comes from the root 
quadruple [— 2, 3, 6, 7 == ]. Isn't "she" lovely! 




Figure 11: Apollonian packing [—2, 3, 6, 7] contains "rectangle" [—2, 3, 10, 15]. 



9 Reverse Trends. 

We have seen how to get many circle configurations out of a single Pythag- 
orean triple. But now we glance at Fig. 12 to discover that one tangent circle 
array can be used to make more triples! This is a simplified version of Fig. 8 
showing the centers C, B, A, D, Z of five circles, each marked with a curvature 
(1,2,3,6,23). 

By various techniques, one can (a) find the center Z of the (23)-circle, 
(b) join Z as shown to the vertices ABC, and (c) drop perpendiculars to 
the sides AB, BC, AC. This divides the triangle ABC into six smaller 
right triangles, two of them congruent (a). In another article [4] we plan to 
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Figure 12: Dissection into rational right triangles. 



analyze the general case which this diagram points to, and show that all the 
small triangles are rational right triangles. Specifically, in this example, the 
triangles a, /3, 7, 5, e are similar to right triangles 

[21,20,29], [5,12,13], [7,24,25], [33,56,65], [117,44,125]. 

The enterprising reader is invited to find appropriate multiples of these 
triangles and fit them together like a jigsaw to make a multiple of [3, 4, 5]. 

10 Table of integer DCE solutions 

The table contains only sorted sequences of form [—a, b, c, d == ] in which d 
is its own reflection (i.e.— a + b, +c = d) or of form [—a, b, c, d ><e ] in which 
d < e and d reflects to e. The first category always gives IAPs that are 
bilaterally symmetric. All entries are "root quadruples" as in [10]. All are 
all obtained by reducing quadruples obtained from PPT's. 

All solutions from (24) are recognizable in that —a + b = d — c = m? 
and - = — . This includes all in the first category. When one inspects a 
diagram such as Fig.ll, the rectangle characteristic of a PT is spotted by 
placing horizontal lines through the centers of the (-2) and (3) circles, joining 
also the (10)s and the (15)s: [—2,3,10,15]. Equivalent numerical tests are 
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Table 1. Same root quads —a, 1j , c, d 



d =: d reflects to d 



_1 


2 


2 


3 




_ 2 


3 


g 






_ 3 


4 


12 


j_3 




_jj 




2 


2 1 




_5 


h 


30 


31 














-6 




42 


43 




_ 7 


g 


5 5 


5 7 




_ g 


y 


7 2 


7 3 




_9 


10 


90 


91 




- 10 


14 


3 5 


3 9 




- 10 




110 


111 




- 11 


12 


13 2 


13 3 




-12 


21 


28 


37 




- 14 


18 


6 3 


6 7 




- 15 


24 


4 


4 9 




- IS 


2 2 


9 9 


10 3 




-20 


36 


45 


61 














-22 


26 


143 


147 




- 24 


3 ^ 




9 7 




- 2 6 


3 


19 5 


19 9 




- 2 S 


44 


7 7 


9 3 




-30 


39 


130 


139 




- 3 


5 5 


5 5 


9 1 




- 3 


34 


2 5 5 


25 9 




- 3 3 


4 2 


154 


16 3 




-35 


50 


84 


109 














-39 


48 


208 


217 




-4 


5 5 


104 


12 9 




-4 2 


5 1 


238 


24 7 




-44 


6 


16 5 


18 1 














-48, 


57 , 


304, 


313 




-51, 


60, 


340, 


349 




"52, 


68, 


221, 


237 




-55, 


80, 


176, 


201 




-57 , 


66, 


418 , 


427 




-60, 


85, 


204, 


229 




-60, 




285, 


301 




-60 , 


59, 


450 , 


459 




-63, 


112, 


144, 


193 




-65, 


90, 


234, 


259 




-66, 


102, 


187, 


223 




-68, 


84, 


357, 


373 




-69, 


78, 


598, 


607 




-70, 


119, 


170, 


219 




-70, 


95, 


266, 


291 




-75, 


34, 


700, 


709 




-76, 


92, 


437, 


453 




-77 , 


12 6, 


19 8 , 


24 7 




-78, 


114, 


247, 


283 




-84, 


133, 


228, 


277 




-85, 


110, 


3 74 , 


399 




-88, 


152, 


209, 


273 




-90, 


115, 


414, 


439 





d X e: d reflects to e 



_3 


5 


g 


8 




12 


_ 5 


'J 


18 


18 




j 2 


- 6 


14 


11 


15 




^ I. 


_ 'J 




3 2 


3 2 




~'< 6 


_9 




50 


50 




54 














-12 


21 


29 


32 




44 


— 14 


2 2 


3 9 


4 3 




5 1 


- 15 




2 8 


4 




5 2 


-15 


24 


41 


44 




56 


— IS 


2 6 


5 9 


R 




7 1 


- 2 


4 5 


3 7 


5 2 




7 2 


- 2 1 


3 9 


4 6 


R ; : ; 




7 


-21 


30 


71 


74 




86 


— 2 2 


3 


8 3 


,:, n 




9 5 


- 24 




8 9 


9 2 




1 04 


— 2 6 


34 


111 


115 




123 


-28 


60 


53 


77 




93 














-30 


39 


131 


134 




146 


- 3 3 


5 1 


94 


10 6 




1 18 


- 3 3 


4 2 


15 5 


15 8 




170 


- 3 5 


84 


64 






14 1 


-35 


60 


88 


93 




133 


— 3 6 


6 8 


7 7 


10 1 




1 17 


- 3 9 


5 7 


124 


13 6 




14 8 


—4 


6 5 


10 8 


113 




153 


-42 , 


91, 


79, 


114 


X 


142 


-44 , 


76 , 


109 , 


117 




165 


-45, 


70, 


127, 


142 


I:-:! 


162 


-45, 


95, 


86, 


126 


X 


146 


-51, 


69, 


196, 


208 




220 


"52, 


84, 


141, 


149 




197 


"55, 


80, 


177, 


192 




212 


-56, 


105, 


121, 


156 




184 


-60, 


92, 


173, 


197 


:=-=: 


213 


-60, 


85, 


208, 


213 


I:-:! 


253 


-63, 


144, 


116, 


161 


I:-:! 


233 


-63 , 


112 , 


148 , 


16 9 




2 25 


-65, 


115, 


154, 


174 


I:-:! 


2 34 


-65, 


90, 


238, 


243 


X 


283 


-70 , 


119 , 


179 , 


18 6 


X 


2 70 


-"70, 


95, 


267, 


282 




302 


-72, 


153, 


137, 


200 




236 


-76, 


108, 


261, 


269 




317 


-77, 


175, 


142, 


198 




282 


-77, 


126, 


207, 


214 




298 


-84, 


133, 


232, 


253 




309 


-84, 


116, 


309, 


317 




365 


-85, 


135, 


234, 


254 




314 


-88, 


209, 


161, 


216 




348 



2(3 



scarce. Certain patterns compel attention. We see the symmetric [—m, m + 
2, n, n >< n+4 ] already discussed. Similarly, we see the pattern [— m,m + 
8, n, n + 4] and solve it: 

m = 2(2Jfe-l), n = 2A; 2 + 2A; - 1. 

One more: [— m, m + 9, n, n + 3] is solved by m = 3k, n — k2 + 3k + 1. We 
have constructed over a dozen such one-parameter families. An added bonus 
is that these formulas have bona fide instances not in the table! Clearly a 
more encompassing theory is to be sought! 
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